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DEVELOPMENT  OF  METHODS  FOR  CALCULATION  OF  THE  THERMAL 
REGIME  OF  ROCK-FILL  DAMS  IN  REGIONS  OF  SEVERE  CLIMATE 

N.A.  Mukhetdinov 

THESIS  ABSTRACT 

The  most  economical  type  of  retaining  structure  in  the  regions  of  the  far  north  is  a dam  made  of 
local  construction  materials.  This  is  true  because  of  the  availability  of  sufficient  quantities  of  local 
construction  materials,  the  possibility  of  using  highly  productive  mechanisms  for  processing  and 
placement  of  these  materials  in  the  structure  with  year-round  continuation  of  construction  opera- 
tions. 

A number  of  water-retaining  structures  have  already  been  built  of  local  construction  materials 
in  Siberia  and  the  far  north,  the  largest  of  these  being  the  dams  of  the  Vilyuyskaya  and  Khantayskaya 
Hydroelectric  Power  Plants.  The  experience  gained  in  the  operation  of  these  dams  has  shown  that 
complex  processes  of  heat  and  mass  transfer  occur  in  them  constantly,  and  determine,  in  the  final 
analysis,  the  formation  of  the  thermal  regime  of  the  structure.  The  temperature  and  moisture  regime 
of  dams  erected  in  regions  with  severe  climatic  conditions  goes  far  to  determine  the  selection  of  the 
type  of  dam  used  and  its  structural  peculiarities. 

In  spite  of  the  successes  which  have  been  achieved  in  the  area  of  planning,  construction  and 
operation  of  dams  of  local  materials,  erected  in  regions  with  severe  climate,  as  was  noted  in  the  reso- 
lutions of  the  coordination  conference  on  water  engineering  construction,  a number  of  problems 
remain  which  require  further  study  before  today’s  requirements  for  planning  and  construction  can 
be  met. 

One  of  these  problems  is  the  development  of  methods  of  calculation  of  the  temperature  and 
moisture  regime  of  rock-fill  dams  considering  the  convection  of  air  and  water  vapor  through  the  pores 
in  the  fill  with  complex  boundary  conditions;  this  is  the  primary  content  of  the  present  dissertation. 

The  work  presents  the  theoretical  and  experimental  (laboratory  and  natural)  investigations  of  the 
author,  revealing  the  processes  of  formation  of  the  thermal  regime  of  individual  structural  elements 
of  a rock-fill  dam  (downstream  prism,  foundation,  central  impervious  core  and  inclined  impervious 
core). 

The  work  was  performed  at  the  Siberian  Affiliate  of  the  All-Union  Scientific  Research  Institute 
of  Hydraulic  Engineering  imeni  B.  Ye.  Bedeneyev  in  1966-1971. 

Chapter  1 studies  the  development  and  current  status  of  theoretical  and  experimental  studies  on 
free  convection  in  porous  media,  the  determination  of  their  heat  and  moisture  exchange  properties. 

A rock-fill  dam  erected  under  difficult  geological,  geocyrological  conditions  of  the  base  is  an  open 
thermodynamic  system.  The  material  and  energetic  interactions  of  the  dam  with  its  environment 
disrupt  its  thermal  equilibrium,  leading  to  the  appearance  of  flows  of  heat,  water,  air  and  vapor  in  the 
dam.  The  intensity  of  the  process  of  heat  and  material  transfer  is  determined,  in  addition  to  the  ex- 
ternal factors,  by  the  physical,  mechanical  and  thermal-physical  properties  of  the  materials  of  the 
foundation  and  the  elements  of  the  dam  (specific  and  volumetric  heat  capacity,  coefficients  of 


conductive,  radiant  and  convective  heat  transfer,  aerodynamic  and  hydrodynamic  resistance  of  the 
fill  material,  temperature  and  moisture  conductivity  coefficients).  There  are  comparatively  numer- 
ous experimental  and  theoretical  works  concerning  all  of  these  problems,  as  well  as  free  air  con- 
vection. 

Analysis  of  the  influence  of  free  air  convection  on  the  thermal  regime  of  a medium  with  large 
pores  forms  the  subject  of  a work  by  M.E.  Aerov  and  N.N.  Umnik.  They  suggested  a specific 
analytical  dependence  for  determination  of  the  effective  (adjusted)  heat  and  temperature  conduc- 
tivity coefficients,  considering  the  intensification  of  heat  transfer. 

V.G.  Gol’dman,  G.V.  Porkhayev,  R.M.  Sarkisyan,  N.V.  Ukhova  and  N.S.  Ivanov  analyzed  free 
convection  in  the  general  process  of  heat  exchange  in  calculating  the  thermal  regime  of  thawing 
permafrost.  However,  the  calculation  formulas  which  they  produced  are  correct  only  for 

(Pages  4 and  5 missing  — tr.] 
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Analysis  of  works  on  the  study  of  free  convection  shows  that  the  general  equations  of  heat  and 
mass  transfer  of  free  [illegible  — tr.]  of  the  air  in  materials  with  large  pores  have  not  yet  been  formu- 
lated, while  the  heat  and  mass  conducting  properties  have  been  studied  basically  only  for  porous 
media  with  small  particle  diameters. 

Further  studies,  both  theoretical  and  experimental,  are  necessary  to  determine  the  regularities 
and  nature  of  changes  in  the  temperature-moisture  regime  of  large-pore  media. 

In  Chapter  2,  the  basic  equations  of  heat  and  mass  transfer  are  formulated  relative  to  the  down- 
stream supporting  prism  of  rock-fill  dams,  and  methods  of  their  numerical  solution  are  analyzed. 

In  the  thermodynamic  respect,  the  downstream  prism  of  a dam  was  analyzed  as  a heterogeneous 
multicomponent  system,  consisting  of  the  mineral  skeleton,  water,  water  vapor  and  pore  gases.  The 
transfer  of  heat  within  this  system  is  by  heat  conductivity,  radiation  and  convection,  the  transfer  of 
matter  — by  molecular  and  molar  movement  of  vapor  and  liquid. 

The  heat  and  mass  transfer  equations  were  produced  on  the  basis  of  the  law  of  conservation  of 
energy  and  matter,  the  law  of  changing  momentum,  under  the  condition  that  the  transfer  coefficients 
can  be  determined  experimentally  with  the  corresponding  potentials. 

The  vapor  space  of  natural  large-pore  media,  particularly  rock  fill,  has  a complex  and  disordered 
structure. 

The  introduction  of  statistical  characteristics  of  porosity,  filtration  rate,  permeability  and  tempera- 
ture allowed  the  flows  to  be  homogenized  within  a small  volume  and  the  concepts  of  speed,  pressure, 
permeability,  temperature  and  mass  at  a given  point  to  be  introduced  to  the  theory  of  heat  and  mass 
transfer.  Therefore,  the  equations  are  correct  for  a fictitious  medium  with  physical  characteristics 
corresponding  to  the  mean  integral  parameters  of  rather  large  volumes  of  actual  fill.  The  boundary 
conditions  of  the  equations  of  heat  and  moisture  transfer  are  strictly  determined  only  if,  in  addition 
to  the  equations  produced,  we  also  analyze  equations  describing  the  speed,  moisture  and  tempera- 
ture fields  of  air  outside  the  dam  considering  the  varying  meteorological  conditions. 

Analysis  of  the  results  of  observations  of  the  dam  of  the  Vilyuyskiy  Hydroelectric  Power  Plant 
over  a number  of  years  has  shown  that  the  boundary  condition  at  the  surface  of  the  downstream 
prism  for  the  equation  defining  the  velocity  field  of  the  air  can  be  determined  with  accuracy  suf- 
ficient for  practice  by  setting  the  normal  derivative  of  the  flow  function  equal  to  0,  the  temperature 
of  the  vapor-air  mixture  at  the  entry  to  the  downstream  prism  equal  to  the  external  air  tempera- 
ture; moisture  transfer  has  no  significant  influence  on  the  formation  of  the  thermal  regime  of  the 
rock-fill  downstream  prism. 

The  system  of  equations  for  determination  of  the  temperature  regime  of  the  downstream  prism 
'of  a rock-fill  dam  can  be  reduced  to  the  form: 
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The  boundary  conditions  are: 


1 .  The  temperature  at  the  permeable  surface  of  the  downstream  prism 


->»,)♦  if-*. 
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where  x = x. 


2.  The  temperature  at  the  interfaces  between  the  downstream  prism  and  the  other  elements  of 
the  dam  (central  core,  inclined  core,  foundation) 

Xk  = Xint  -jjp  J fint  * 1 • W 


3.  The  temperature  of  the  moving  air  at  the  boundaries  of  the  downstream  prism 


a)  0{x,  y,  t)  = 0(x'sur,y'sur,  t)  = fwt 

where  x = x'wr 
y =>'sur  I 

b)  0(x,  y,  t)  = e[xint,ylnt,  r)  = t(xinvyinv  r), 


where  x = xjnt 

y=y>nf 

4.  The  flow  line  function 


a)  4>-0, 

where  0<x<fl,y  = 0and0<>'</y,x  = 0 

where  x = x^r  \y  -yw ; 
c)  ^ a const 


in  the  air-impervious  sectors  of  the  downstream  prism. 
The  initial  conditions  are: 


tb.y,  r)=t(x,y,0); 


r)  =6(x,y,0); 

'hiy.y.r)  = 4>{x,y,  0)  (9) 

where  t,  6 are  the  mean  integral  temperature  of  the  rock  fill  and  of  the  air,  °C; 

fsUri  ®Sur are  the  same,  at  the  pervious  surface  of  the  downstream  prism,  °C; 
rext  is  the  temperature  of  the  external  air,  °C; 

^ is  the  flow  function,  M2/hr; 

dv  is  the  volumetric  heat  exchange  coefficient,  kcal/M3  hr  °C; 

0 is  the  coefficient  of  volumetric  expansion  of  the  air,  1 1° C; 
p is  the  density  of  the  air,  kg/M3 ; 
p0  is  the  density  of  the  air  at  0°C,  kg/M3 ; 
y is  the  volumetric  weight  of  the  air,  kg/M3; 
m is  the  porosity  of  the  fill; 

Cvoi,  ^voi 's  t^,e  volumetric  heat  capacity  of  the  fill  and  air,  kcal/M3  °C; 
q is  the  quantity  of  heat  carried  into  the  dam  by  heat  conductivity,  kcal/M3  hr; 

Xk  (dr/d/?)  is  the  heat  flux  from  the  downstream  prism  to  the  line  of  contact  with  elements 
of  the  dam,  kcal/M2  hr; 

Xint  (dtim/dn)  is  the  same,  from  the  central  core  (inclined  core),  and  base  of  the  dam, 
kcal/M2  hr; 

Xk  is  the  coefficient  of  transverse  convective  heat  conductivity,  kcal/M  hr  °C; 

Xj„t  is  the  same,  for  the  base,  central  core  (inclined  core),  kcal/M  hr  °C; 
n is  a perpendicular  to  the  surface,  M; 

xinvym  are  the  coordinates  of  the  interface  of  the  downstream  prism  with  the  base,  central 
core  (or  inclined  core),  M; 

Xjur^jur are  the  coordinates  of  the  surface  of  the  downstream  prism,  M; 

*’sur>y'sur  are  the  coordinates  of  the  surface  of  a fictitious  layer,  M; 

B is  the  width  of  the  downstream  prism  at  the  base,  M; 

H is  the  height  of  the  downstream  prism,  M. 

Equation  system  (1 ) was  used  to  calculate  the  thermal  regime  of  an  isotropic  and  anisotropic 
rock-fill  downstream  prism. 

In  the  isotropic  prism,  considering  the  nonlinearity  of  filtration  of  air  with  free  convective  motion, 
the  coefficients  of  equations  (1 ) and  (3)  were  taken  as: 

CVO|  <t>y)  = Cvoi  ■ con*1: 
m(x,  y)  =m  - const; 

X(x,y)  = Xef  = const. 

jx/p*,  at  Re  = Rea 

Kx=Ky  - f{f>w)  V (10) 

1 1/p  [{ulk)  * B0pw),  at  Re  > Reer 
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where  Xef  is  the  coefficient  of  heat  conductivity  of  the  fill,  kcal/M  hr  °C; 

H is  the  dynamic  viscosity  factor  of  the  air,  kg  hr/M2 ; 

K is  the  coefficient  of  permeability  of  the  fill,  M2; 

B0  is  the  turbulence  parameter,  kg/M4; 

Re  is  the  Reynolds  number; 

Re^  is  the  critical  Reynolds  number,  defining  the  transition  of  the  air  flow  from  one  mode 
to  the  other. 

2.  In  an  anisotropic  downstream  prism,  the  coefficients  of  equation  (3)  are 


/0f=/Cy= A ; 

m =m(y) ; 

^voi - ^voi  (*f  y)- 


Equation  system  (1)  includes  nonlinear  equations  of  parabolic  (1),  hyperbolic  (2)  and  elliptical 
(3)  types.  Solution  of  the  system  of  nonlinear  equations  in  analytic  form  is  possible  only  in 
exceptional  cases.  The  universal  method  of  approximate  solution  of  differential  equations  is  the 
method  of  finite  differences,  which  was  used  in  performing  our  calculations.  The  author  has  de- 
veloped a special  method  of  numerical  solution  for  the  quasilinear  elliptical  equation  describing 
the  motion  of  air  considering  the  nonlinearity  of  filtration. 

The  equation  of  nonlinear  filtration  (3),  considering  (1 0),  is  reduced  to  the  form: 

where  pQ,  p0, 0O  are  the  density,  viscosity  and  coefficient  of  volumetric  expansion  of  the  air  cor- 
responding to  the  average  temperature  of  air  moving  through  the  pores  over  the  profile  of  the 
downstream  prism. 

Let  us  introduce  the  function 


4>{x,y,  r) 


Po0o9rr> 


Equation  (1 2)  is  transformed  to 

92< t + d20  _ _ d6_ 

3x2  by2 

This  equation  should  be  solved  with  the  boundary  and  initial  conditions 
0 = 0,  y*0, 0<x<8;x30,  0<yc/f; 


= 0,  where  x = x^r,y  -yw\ 


<t>(x,y,  r)=0(x,y,  0). 


The  numerical  solution  of  equation  (13),  considering  boundary  conditions  (14),  is  usually  pro- 
duced by  finite-difference  methods  using  a plan  of  traveling  calculation  and  simple  iteration.  The 
solution  of  this  equation  was  considered  to  be  found  if  the  difference  in  two  successive  approxima- 
tions satisfied  the  condition 

|i«£|i-i*S,j||  ««.  (15> 

where  S is  the  number  of  iterations; 

e is  a predetermined  comparison  constant. 

Based  on  the  known  values  of  function  <t>(x,  y,  r)  at  discrete  points  in  the  area  in  question,  die 
following  algebraic  equations  were  produced  for  determination  of  the  flow  function: 


* T i * [ (2<'lcL' ' ~ 


2-  *1  *.j  * £ 1 2%^  - - vhj  * - . (”) 

where  v0  is  the  coefficient  of  kinematic  viscosity  of  the  air  corresponding  to  the  average  tempera- 
ture of  air  moving  through  the  pores  across  the  profile  of  the  downstream  prism,  M2/hr; 

c=/-i_+  JL-V 

l A*2  Ay2/' 

V A*2  £y2  / 

a = / ^2=14-  + . (18) 

\ Ax2  Ay2  / 

Arsh  is  the  reciprocal  value  of  the  main  hyperbolic  sine. 

The  stability  and  convergence  of  the  finite-difference  plans  for  equation  system  (1 ) are  defined 
from  the  conditions: 

1 _ *.Xn-lj  + + dyt„\  —AT..  I > 0 

* ' ^vol.n.i  J 

where  X„  ; Xn  j_, ; Xn+l  j;  Xn_,  j is  the  heat  conductivity  coefficient  of  points  with  the  correspond- 
ing indices,  kcal/M  hr  °C. 
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Chapter  3 presents  the  results  of  experimental  studies  of  physical-mechanical  and  heat-physical 
properties  of  a rock  fill. 

Processing  of  the  results  of  laboratory  experiments  and  field  studies  allowed  the  production  of 
empirical  equations  for  determination  of: 

- the  surface  area  (F)  of  an  individual  unit  of  fill  volume 

F (20) 

d ’ 

where  d is  the  calculated  mean  diameter  of  an  individual  rock,  M; 

- the  mean  characteristic  dimension  (/?v)  of  the  fill 


/?„  = 


v~  13.5  (1  -m) 


- the  aerodynamic  resistance  (fp)  of  the  large  pore  fill 


(21) 


(22) 


- the  coefficient  of  permeability  (fo)  (where  d > 10  cm) 


K=  d2-mi-  - , (23) 

3100  (1-m)2 

- the  turbulence  parameter  ( B0 ) 

B0  = 0.83  . (24) 

dm3 


The  heat  conductivity  coefficient  of  the  fill  (Xef)  is  satisfactorily  described  by  the  formula  of 
V.l.  Odelevskiy. 

The  temperature  conductivity  coefficient  ( a 0)  is  defined  as  the  derivative  of  a quantity: 


a 


-h* 


(25) 


l-i 
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The  volumetric  heat  capacity  (Cvol)  of  the  fill  was  found  using  the  formula 

Cyoi  = (1  -m)  Ct p + mC^  , (26) 

where  Cv  is  the  specific  heat  capacity  of  the  material  of  the  rock  fill,  kcal/kg  °C; 
yck  is  the  volumetric  weight  of  the  material  of  the  rock  fill,  kg/M3. 

Equations  (20-26)  allow  us  to  determine  all  of  the  parameters  of  the  two-phase  rock  fill  which 
we  will  need. 

In  Chapter  4,  results  are  presented  from  calculation  of  the  influence  of  nonlinear  filtration  of 
air  on  the  formation  of  the  thermal  regime  of  the  rock-fill  downstream  prism  and  the  calculated 


values  of  temperature  are  compared  with  the  data  of  field  observations  of  the  rock-fill  dam  of  the 
Vilyuyskaya  Hydroelectric  Power  Plant. 

In  determining  the  thermal  regime  of  the  dam  of  the  Vilyuyskaya  Hydroelectric  Power  Plant, 
two  structural  models  of  the  downstream  prism  were  analyzed: 

a)  Anisotropic,  when  the  mean  calculated  diameter  of  individual  rocks  and  the  porosity  of  the 
fill  in  the  vertical  plane  were  taken  as  functions  of  coordinates,  but  they  were  considered  constant 
in  any  horizontal  plane; 

b)  Isotropic,  with  porosity  and  diameters  of  individual  rocks  corresponding  to  the  mean  values 
of  the  anisotropic  fill. 

The  upstream  shell,  foundation,  and  inclined  core  were  looked  upon  as  a rigid  body  with  varying 
heat-physical  characteristics,  in  which  heat  transfer  occurs  only  by  conduction.  The  latent  heat  of 
phase  transitions  during  freezing  (and  thawing)  in  the  inclined  core  and  upstream  shell  were  con- 
sidered by  a method  developed  by  Yu.  Ye.  Gavrish  and  L.N.  Khrustalev.  Filtration  of  air  in  the 
downstream  prism  was  considered  by  the  rule  of  Darcy. 

The  filtration  anisotropy  of  the  downstream  prism  of  the  dam  leads  to  redistribution  of  the  air 
flow,  qualitative  and  quantitative  changes  in  its  thermal  regime,  particularly  in  areas  near  the  slope 
and  foundation  of  the  dam.  In  an  anisotropic  fill,  in  comparison  with  an  isotropic  fill,  the  concen- 
tration of  flow  lines  at  the  foundation  shows  that  the  main  paths  of  movement  of  air  in  the  down- 
stream prism  are  sections  of  the  fill  with  large  diameters  of  individual  rocks  and  high  porosity.  As 
we  move  away  from  the  surface  of  the  slope,  the  influence  of  anisotropy  on  the  formation  of  the 
thermal  regime  of  the  downstream  prism  decreases.  In  the  zone  beneath  the  inclined  core,  the  tem- 
perature fields  produced  by  calculations  for  both  media  are  almost  identical  or  differ  little  from 
each  other,  and  also  agree  closely  with  the  observed  data. 

The  difference  in  temperature  between  the  calculated  and  observed  data  was:  in  the  zone  be- 
neath the  inclined  core  ±2°C  and  in  the  surface  zone  of  the  downstream  prism  +4°C.  The  main 
reason  for  the  difference,  as  observations  showed,  was  the  diurnal  and  weekly  fluctuations  in  tem- 
perature, leading  to  changes  in  the  direction  of  air  flow  within  the  body  of  the  downstream  prism. 

The  predominant  direction  of  natural  movement  of  the  air  flow  in  the  downstream  prism  is  deter- 
mined by  the  external  air  temperature.  In  the  water,  the  flow  enters  at  the  base  and  leaves  through 
a sector  near  the  top  of  the  dam,  while  in  the  summer  the  air  flow  is  in  the  opposite  direction.  In 
accordance  with  this,  in  moderate  or  high  head  dams,  in  the  upper  half  of  the  downstream  prism 
the  temperature  of  the  fill  changes  its  sign  during  the  course  of  the  year,  while  in  the  lower  half  of 
the  downstream  prism  in  areas  with  year-round  negative  temperatures,  the  temperature  remains 
negative  [below  freezing]  constantly.  In  the  fall,  during  a period  of  1 .5-3  months,  we  see  a shift 
in  the  isotherms  of  negative  temperature  toward  the  impervious  elements  of  the  dam. 

The  rate  of  change  of  the  thermal  regime  of  the  downstream  prism  can  be  regulated  by  placement 
of  an  air-impervious  blanket  on  its  slope.  The  primary  effect  of  the  air-impervious  blanket  is  that  it 
excludes  convective  heat  exchange  in  the  downstream  prism  with  the  atmosphere.  The  rate  of  move- 
ment of  the  air  in  the  pores  is  greatly  decreased,  the  mean  temperature  of  the  fill  changes  more 
slowly  than  in  dams  with  free  air  exchange. 

In  Chapter  5,  a study  is  made  of  individual  problems  of  the  thermal  regime  of  the  core  and  founda- 
tion of  a rock-fill  dam.  The  core  of  a rock-fill  dam  experiences  the  influences  of  complex  processes 
of  heat  exchange,  occurring  at  the  surface  of  the  dam,  in  the  supporting  downstream  prism,  and  also 
in  the  upstream  shell,  which  is  in  thermal  interaction  with  the  reservoir.  In  planning  and  construction 
of  rock-fill  dams,  it  is  important  to  establish  the  general  trend  of  formation  of  the  thermal  regime  of 
the  core  and  foundation  of  the  dam,  if  possible  considering  all  of  the  factors  mentioned. 

The  formation  of  low  negative  temperatures  in  the  fill  of  the  downstream  prism  and  their  mainte- 
nance throughout  the  year  creates  conditions  for  freezing  of  the  foundation  of  the  rock-fill  dam. 
Depending  on  the  dimensions  of  the  downstream  prism,  the  “aureole"  of  freezing  of  the  nonfiltering 
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foundation  of  the  dam  may  become  significant.  It  can  be  determined  quantitatively  by  the  method 
of  calculation  of  thawing  beneath  industrial  and  civil  structures. 

Under  more  complex  conditions  of  mass  transfer,  a zone  of  negative  temperature  is  formed  in 
filtering  foundations,  and  depends  not  only  on  the  size  of  the  downstream  prism  and  temperature 
at  the  base  (rb),  but  also  on  the  filtration  factor  (/Cf),  the  heat  and  temperature  conductivity  of  the 
soil  of  the  base  when  thawed  and  frozen  (Xt,  ak  Xm ),  the  head  (H)  and  temperature  of  the  bottom 
layers  of  water  in  the  reservoir  (twU).  However,  the  formation  of  an  “aureole”  of  freezing  in  the 
foundations  of  rock-fill  darm  creates  the  effect  of  a cut-off  in  relationship  to  the  filtration  flow. 

This  decreases  the  filtration  flow  passing  through  the  foundation. 

The  maximum  freezing  depth  of  the  foundation  of  the  dam  with  filtration  heat  transfer  is  dis- 
placed in  the  downstream  direction  and  can  be  determined  from  the  equation 

Smax  -f(v  \ 

81  ~ \ VwatVM  / 

where  £max  is  the  maximum  depth  of  freezing  of  the  foundation  beneath  the  downstream  prism,  M; 
2,  is  half  the  length  of  the  horizontal  subterranean  contour  of  the  dam,  M. 

Studies  have  shown  that  the  temperature  of  the  filtering  flow  in  the  foundation  at  the  lower  face 
of  the  core  differs  very  little  and  is  always  lower  than  the  temperature  of  the  water  in  the  reservoir. 

Reverse  filters,  in  which  heat  is  propagated  only  by  conduction,  significantly  decrease  the  cooling 
influence  of  the  downstream  prism  on  the  impervious  elements. 

The  minimum  thickness  of  a reverse  filter  providing  for  exit  of  filtration  water  at  0 head  can  be 
determined  by  the  following  approximate  formula: 


W] 


^^^wat^wat 


where  2 is  the  thickness  of  the  reverse  filter,  M; 

is  the  effective  heat  conductivity  coefficient  of  the  material  of  the  reverse  filter,  kcal/M 

hr  °C; 

in  is  the  minimum  temperature  of  the  fill  of  the  downstream  prism  at  the  boundary  with 
the  reverse  filter,  °C; 

X3  is  the  heat  conductivity  factor  of  the  first  layer  of  the  filter  in  the  thawed  state,  kcal/M 

hr  °C; 

G is  the  flow  rate  of  water  at  the  line  of  the  surface  of  the  core,  kg/M2  hr; 
fWJI  is  the  temperature  of  the  water  in  the  reservoir,  °C; 

Cwat  is  the  specific  heat  capacity  of  the  water,  kcal  kg/°C; 

h is  the  distance  from  the  surface  of  the  water  filtering  through  the  downstream  face  of  the 
core  to  the  0 isotherm,  M. 


Under  the  cooling  influence  of  the  downstream  prism,  the  core  is  frozen  (filtration  is  stopped) 
in  dams  in  which  the  flow  rate  of  water  through  the  core  is  less  than  that  calculated  by  the  formula: 


(^-wat  ^wat  + 1<t>^ 

where  is  the  latent  heat  of  phase  conversion,  keal/kg. 
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CONCLUSIONS 

1 . Physical  equations  are  produced  for  the  heat  and  mass  transfer  in  the  downstream  prism  of  a 
rock-fill  dam,  and  conditions  of  unambiguity  of  these  equations  are  established. 

2.  Based  on  actual  observations,  the  possibility  is  established  of  simplification  of  the  boundary 
conditions  of  equations  describing  the  temperature-moisture  regime  of  the  downstream  prism. 

3.  Quantitative  relationships  are  established  for  the  change  in  surface  area  of  individual  rocks  in 
a unit  of  fill  volume,  the  mean  characteristic  diameter,  aerodynamic  resistance  and  heat-physical 
characteristics  of  the  rock  fill. 

4.  A method  is  developed  for  calculation  of  the  temperature  fields  of  rock-fill  dams  considering: 

- The  thermal  and  filtration  anisotropy  of  the  downstream  prism; 

- The  conditions  of  heat  exchange  of  a rock-fill  downstream  prism  with  the  surrounding  environ- 
ment, variable  in  time  and  space; 

- Heat  transfer  of  the  filtering  flow  of  water  in  the  foundation  and  the  core  of  the  dam; 

- The  engineering  geological,  geocryological  conditions  of  the  foundation  and  the  geometric 
forms  of  dam  elements; 

- Nonlinearity  of  filtration  with  free  convective  motion  of  air  in  large-pore  media. 

5.  An  algorithm  has  been  written  and  calculations  performed  on  a computer  to  determine  the 
thermal  regime  of  the  rock-fill  dam  of  the  Vilyuyskaya  Hydroelectric  Power  Plant  for  anisotropic 
and  isotropic  models  of  the  downstream  prism. 

The  results  of  the  calculation  and  the  data  of  actual  observations  agree  closely. 

6.  The  anisotropy  of  the  downstream  prism,  resulting  from  the  technology  of  work  involved  in 
placing  the  rock  fill  in  layers,  leads,  in  comparison  to  an  isotropic  prism,  to  qualitative  and  quanti- 
tative redistribution  of  the  air  flow  and  temperature  in  the  prism  only  in  areas  near  the  downstream 
slope. 

7.  Intensive  change  in  the  thermal  regime  of  the  downstream  prism  of  rock-fill  dams  occurs  during 
the  first  two  months  after  the  outside  air  temperature  falls  below  the  freezing  point. 

8.  In  rock-fill  dams  erected  in  severe  climatic  regions  with  free  air  exchange  between  the  down- 
stream prism  and  the  atmosphere,  two  characteristic  thermal  regime  zones  are  formed: 

- The  upper  portion  of  the  downstream  prism,  where  the  temperature  rises  above  and  falls  below 
freezing  during  the  year; 

- The  lower  portion  of  the  downstream  prism,  where  the  temperature  remains  below  freezing 
throughout  the  year. 

9.  An  air-impervious  blanket  on  the  surface  of  the  downstream  prism  has  a significant  influence 
on  the  overall  thermal  regime  of  rock-fill  dams.  The  absence  of  convective  heat  exchange  between 
the  downstream  prism  and  the  outside  air  during  the  summer  results  in  cooling  and,  during  the 
winter  - in  warming  of  the  overall  temperature  of  the  rock-fill  dam. 

10.  Calculations  have  shown  that: 

- The  natural  thermal  regime  of  a rock-fill  downstream  prism,  even  if  filtration  is  eliminated,  does 
not  provide  the  freezing  of  all  impervious  elements  of  a frozen  type  dam  required  for  safe  use; 

- Elimination  of  periodic  freezing  of  the  core  in  dams  of  this  type  and  free  passage  of  water 
filtering  through  the  core  can  be  assured  by  proper  selection  of  the  necessary  thickness  of  reverse 
filters. 

1 1 . Studies  of  the  freezing  of  the  filtering  foundation  of  a rock-fill  dam  have  established  that: 

- In  dams  with  moderate  or  high  head,  the  temperature  of  the  water  filtering  into  the  foundation 
beneath  the  central  core  (or  inclined  core)  is  close  to  the  temperature  of  the  water  in  the  reservoir; 

- The  maximum  freezing  depth  of  the  foundation  of  the  dam  is  displaced  from  the  downstream 
slope  by  a distance  of  2/3  of  the  half  length  of  its  underground  contour; 

- The  freezing  of  the  bed  beneath  the  downstream  prism  of  a rock-fill  dam  significantly  decreases 
or  even  eliminates  filtration  in  the  foundation  of  the  structure. 
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